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Abstract: We calculate the quark coefficient function T'^{x,^) that enters the factorized 
amplitude for deeply virtual Compton scattering (DVCS) at all order in a soft and collinear 
gluon approximation, focusing on the leading double logarithmic behavior in (xziz^), where 
X lb ^ is the light cone momentum fraction of the incoming/outgoing quarks. We show 
that the dominant part of the known one loop result can be understood in an axial gauge 
as the result of a semi-eikonal approximation to the box diagram. We then derive an all 
order result for the leading contribution of the ladder diagrams and deduce a resummation 
formula valid in the vicinity of the boundaries of the regions defining the energy flows of 
the incoming/outcoming quarks, i.e. x = ib^. The resummed series results in a simple 
closed expression. 
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1. Introduction 

Since a decade, there has been much progress in the understanding of the three-dimensional 
content of the hadron, both from the theory and the experimental sides. Experimentally, 
this relies on several new electron facilities combining high luminosity and advanced detec- 
tors which allow for measuring with an impressive precision exclusive processes, including 
deep virtual Compton scattering (DVCS) and meson production. This lead to the first 
studies of non-perturbative non forward parton distributions, now called generalized par- 
ton distributions (GPDs), first in the fixed target experiment HERMES [1-3], and then 
at HI and ZEUS, using the dominance of the DVCS contribution at small xsj [4-6]. Al- 
most simultaneously, the DVCS contribution was measured at JLAB, at CLAS [7,8] and 
at Hall A [9]. Prom the theory side, the interest for hard exclusive processes started with 
the Leipzig group [10]. Several studies^ then set the basis of a consistent framework, 
called collinear factorization, to separate the short distance dominated partonic subpro- 
cesses and long distance hadronic matrix elements, at leading and next-to-leading order 
for DVCS [15-20] and for hard electroproduction of mesons [21-24] and their timelike 
crossed versions, namely exclusive lepton pair production in photon or meson collisions 
with protons [25-27]. The future JLab-12 GeV and COMPASS-H program will provide 
soon bunches of data, giving a hope to get access to GPDs with a high degree of precision. 
There are indeed now intense activities [28-31] to move from discovery era to precision 
physics. 

In order to extract the GPDs, a precise theoretical framework should be available, 
which should go beyond a pure leading logarithmic treatment both for evolution equations 
and for coefficient functions, with the expected increase of precision of future data. The 
aim of this paper is to study in detail the emergence of the leading contributions near 
the points x = ±^ and to derive a resummed formula for the coefficient function of DVCS, 
which could have a major phenomenological impact in future precise studies. A brief report 
on this result has been presented elsewhere [32]. 

The amplitude for the DVCS process 

7^*Hq)N{p)^^{q')N\p'), (1.1) 

with a large virtuality = —Q^, factorizes at the leading twist 2 level in terms of pertur- 
batively calculable coefficient functions C(x,^, a^) and GPDs F{x,^,t), where the scaling 
variable in the generalized Bjorken limit is the skewness ^ defined as 

s = 91 

{p + p')-{q + q') 

Hereafter, we only consider quark exchange. After proper renormalization, this quark con- 
tribution to the symmetric part of the factorized Compton scattering amplitude illustrated 

^For reviews, see Refs. [11-14] 
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Figure 1: Factorization of the DVCS amplitude in the hard regime. The crossed-blob denote a 
set of r matrices. In this paper F, = i^, . In the above (hard) part, cahed coefficient function, the 
lines entering and exiting the crossed blob carry spinor and color indices but do not propagate any 
momentum. The corresponding momenta are on-shell. 



in Fig. 1 reads 
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where the quark coefficient function T"? read [33] : 
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The first (resp. second) terms in Eqs. (1.5) and (1.7) correspond to the s— channel (resp. 
It— channel) class of diagrams. One goes from the s— channel to the «— channel by the 

interchange of the photon attachments. Since these two contributions are obtained from 
one another by a simple [x -H- —x) interchange, we will restrict in the following mostly to 
the discussion of the former class of diagrams. 
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Eqs. (1.5) and (1.7) show that among the corrections of 0(qs) to the coefficient function 
the terms of order [log^(^±x)]/ {x±^) play an important role in the region of small (^±x), 
i.e. in the vicinity of the boundary between the so-called ERBL and DGLAP domains where 
the evolution equations of GPDs take distinct forms. We here scrutinize these regions and 
demonstrate that they are dominated by soft-coUinear singularities. 

The source of these singularities can be understood in the following way. In our analysis 
we expand any momentum in the Sudakov basis pi, P2, where p2 is the light-cone direction 
of the two incoming and outgoing partons {p^ = p\ = 2pi ■ P2 = s = as 

k = api+pp2 + k±. (1.8) 

In this basis, 

q = Pi-2^P2 and pi = q' . (1.9) 

Now, the Mandelstam variables S and U for the coefficient function illustrated in the upper 
part of Fig. 1 read 

5 = ^Q^ and U = -^Q\ (1.10) 

Although the usual collinear approach is based on a single-scale analysis, where the only 
large scale is provided by Q^, in the special kinematical limit where x — >■ ^ (resp. x — > — ^), 
the Mandelstam variable <S (resp. U) becomes parametrically small with respect to Q^. 
We thus turn to a two scale problem, in a similar way as for the x xbj limit of deep 
inelastic scattering (DIS) on a parton of longitudinal momentum fraction x. In this limit, 
large terms of type [a^ log^(^±x)]/(xib^) should appear, calling for a resummation of these 
threshold singularities, similarly to the resummation of large xbj/x coefficient functions in 
DIS [34,35]. As for DIS, the resummation which we now perform is due to the combination 
of soft and collinear singularities. The main complication with respect to DIS is due to the 
non forward kinematics of DVCS. Our treatment relies on a diagramatic analysis, which 
we now explain. 

We start our analysis by observing that in the same spirit as for evolution equations, 
the extraction of the soft-coUinear singularities which dominate the amplitude in the limit 
X — > ±4 is made easier when using the light-like gauge pi ■ A = 0. We argue that in this 
gauge the amplitude is dominated by ladder-like diagrams, illustrated in Fig. 2. 

We now restrict our study to the case x ^ The dominant kinematics is given by 
a strong ordering both in longitudinal and transverse momenta, according to 

x~^> - |.'r-^| » |x-C + /3i| ~ I/32I > ••• 

••• » |x-e + /3i + /32---- + /3n-i| ~ |/3n|, (1.11) 

l^iil < \kl2\ < • • • < \kin\ < s ~ Q^ (1.12) 

|ai| < ••• < |a;„| < 1, (1.13) 

where and /3j are momentum fractions along the two dominant light cone directions of 
the exchanged gluons and k±i their transverse momenta. This ordering is related to the 
fact that the dominant double logarithmic contribution for each loop arises from the region 
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Figure 2: The ladder diagrams which contribute in the Ught-Uke gauge to the leading a" ln^"(^ — 
x)/{x — ^) terms in the perturbative expansion of the DVCS amplitude. The p2 and _L momentum 
components are indicated. The dashed lines show the dominant momentum flows along the p2 
direction. 



of phase space where both soft and collincar singularities manifest themselves. In the limit 
X ^ ^ , the left fcrmionic line is a hard line, from which the gluons arc emitted in an eikonal 
way, with a collinear ordering. For the right ferniionic line, an eikonal approximation is not 
valid, since the dominant momentum flow along p2 is from the gluon to the fermion. Even 
though this is the case, a collinear approximation can still be applied. This non-symmetric 
treatment of the whole diagram will be referred to as the semi-eikonal approximation. 

The issue related to the ie prescription in Eq. (1.7) is solved by computing the coeffi- 
cient function in the Tinphysical region ^ > 1. After analytical continuation to the physical 
region < ^ < 1, the physical prescription is then obtained through the shift ^ ^ — ie . 

We denote by Kn the contribution of a n-loop ladder to the coefficient function, and 
define^ 

1 / 1 

As mentioned earlier, Eq. (1.14) should be completed by inclusion of the u-channel class 
of diagrams. 

^Note that this normalization is fixed to match with the definition of the coefficient function used in 
Eq. (1.4), while the 5-matrix element is iKn at order n. 
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Figure 3: One-loop self energy diagram. 

The paper is organized as follows. In Sec. 2, we perform a detailed analysis of the one- 
loop diagrams and extract the dominant contribution in [og log^{^-x)]/{x-^). In Sec. 3, we 
analyze the same one- loop contribution but in the spirit of the semi-eikonal approximation 
described above and show that the dominant contribution is indeed identical. In Sec. 4, we 
analyze the two loop-contributions and deduce four guiding rules which are used in Sec. 5, 
in order to demonstrate that only ladder-like diagrams are responsible for [a^log^'^(^ — 
x)]/ {x — ^) contributions, which we then compute and resum. We end up with conclusions 
in Sec. 6. Two appendices give technical details on the analysis of the pole positions 
entering the loop integrals, and on integrals used to extract the dominant contributions. 



2. One- loop analysis based on Ward identities 

In this section we analyze the one-loop diagrams in details without making any approxi- 
mation to understand which diagrams give contribution at order [ag log^(^ — x)]/{x — ^) 
and which give less singular contributions in light-like gauge. We explicitly show that the 
net contribution to [log^(^ — x)]/{x — terms arises from the box-diagram in the case of 
cutting the gluonic line. Moreover, this analysis precisely identifies the part of the phase 
space that is responsible for this contribution. 

2.1 Self energy 

Let us start with the self-energy diagram, illustrated in Fig. 3. The numerator reads 
(Num)s.E. = tr|;^27l [^1 + {x - {i\ + - i)h " ^\ l"" + (^^ - ^^^2] 7±a 

After some algebra, one realizes that the gauge part of the numerator vanishes 



(2.1) 



(Num) 



gauge 



. 



(2.2) 



Since this is the case, the self energy diagram is exactly the same in Feynman and in 
light-like gauge. This diagram is calculated in Feynman gauge in Ref. [33] and it is shown 
that only single log's arise. Hence, the self-energy diagram does not contribute to ^"^^ff^^^-* 
terms. 
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Figure 4: Right-vertex diagram. 
2.2 Right vertex, left vertex and box diagram 

The numerator for the right- vertex diagram shown in Fig. 4 in the hght-hke gauge is written 
as 



(Num)R.v. = tr|;^27'' + {x - Oh] ll + ^ + - Oh] 7^^ [h + - Oh] 7± 

{ k-pi 

After some algebra, one gets 

(Num)R.v.g^. = -8^2 (1 + + X - , 
(Num)R.v. gauge = 8s(/3 + X-0 
Hence, the whole numerator reads 



fc2 "1 



(Num)R.v. =8s(;3 + a;-0^■ 



.(2.3) 

(2.4) 
(2.5) 

(2.6) 



Hereafter, k_\_ (resp. fe) denotes the transverse component of the gluon momentum in 
Minkowski (resp. Euclidean) space and we use = whenever it is needed. Hence, 
the integral for the right-vertex diagram is^ 

s f 1 II 1 

-^R.V. = -2 J dad(5d2k 8g— (^ + x-,^) _ 



[k + {x- Op2)] ^k^[k+p^ + {x- Op2)] ' ■ 

(2.7) 

Similarly, the numerator for the left-vertex diagram in light-like gauge is written as 



(Num)L.v. = tr <j hll [h + {x- Oh] l"" + + (x - OH + (a^ + K)H 

x{,,.- '-^^- + ^-^^- }.(2.8) 



^The angular integration is straightforward for our calculation, which we emphasize through the notation 
djvfc — \k^^^d\k\. This angular integration gives a factor of (27r)" for n-loops, from which 7r" is included 
inside the parenthesis of Eq. (1.14) and 2" is gathered inside /„. 
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Figure 5: The one-loop box diagram. 



After some algebra one gets 



(Num)L.v.ff^. = -8s2(^ + X + 0(1 + a) , 



(Nuin)L.V. gauge = -^{P + X + ^ 



2s2/3(i + a) + 2skj 



(2.9) 
(2.10) 



Hence, the whole numerator reads 



fc2 



(Num)L.v. = 8s-^(^ + x + • 



(2.11) 



One can write the integral for the left-vertex as 



II. 



V, 



-If 



da dp d2k 



1 



8s=r-(P + X + C)—^— TT^ 

. P s{x-01[k+{x + 0P2)] '' [k+Pl + {x-^)p2)] 
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(2.12) 



Let us now calculate the numerator for the box diagram, illustrated in Fig. 5, which is 



(Num)box = tr <j ^2Y [^ + {x- Oh] il [^+h + ix- Oh] l±a [it + {x + Oh] 1^ 

kjiPxv + Kp\^l \ 



9ixv 



k-pi J ■ 



(2.13) 



The Qfj^i, part of the numerator reads 



(Num)boxg,. = tr|[^ + (x+0;^2]7'^MM[^ + (^-0;^2]7l[^+;5i + (^-0;52]7±4 . (2.14) 
Using ^^hlii. — ~'^hi it be written as 



(Num)boxs,. = -2tr []^ + {x + Oh]h[^ + {x - Oh]ll[^ + h + {x - Oh]l ±A ■ (2-15) 



Applying the Ward identity by noting that p2 inside the trace can be put in the form 



p'i = Y^[[k + {x + OP2] -[k + {x- OP2] 



(2.16) 



7 



one gets 



+1 [k+{x- 0P2] '1^1 - (/3 + X + e)| + ^as|. (2.17) 

In this way, we effectively reduce the box diagram to right and left vertex diagrams. The 
gauge part of the numerator is written as 

(Num)boxga«ge = + {x + OM^Mi^ + - OMll + ^ + - Ohh±a 

--^trj [^ + ix + Oh] mh [^ + {x- Oh]il [^ + h + {x- Oh]i^a^ ■ (2-18) 

Since pi = Oj can add or subtract ^2 from ^ when they are appearing next to each 
other inside the trace without spoiling the result, which allows to cancel one of the fermionic 
propagators. The whole numerator for the box diagram is given as 

i2 r 1 



(Num)box = 8[k + {x-OP2_ 
- 8[k + {x + 0P2_ 



( 

2jl 



fci-(^ + x + C)-+^as 



-{l + a){p + x + 



kl-{P + x-Oi^ 



-{l + a){l3 + x-0\. (2.19) 



Rewriting the transverse component of the gluon momentum in Euclidean space, the inte- 
gral for the box diagram is 



'box 



-1/ 



da dp d2k 8 



1 



k^ + iP + x + 07:-(as 



--{l + a){P + x + 



1 



k^ 



J dadpdik (-8)|^ 



[k + {x + OP2] [k+pi + {x- OP2] ' 



(2.20) 



k^ + {P + x-0\ 



+ -{l + a){P + x-0 



k"^ [k+pi + {x- OP2] [k+ [x - OP2] 



2 • 



(2.21) 



Note that the term (2.20) is effectively the same as the left-vertex diagram since the 
fermionic propagator on the outgoing quark line is cancelled. On the other hand, the term 
(2.21) is effectively the same as the right- vertex diagram since the fermionic propagator on 
the incoming quark line is cancelled. This decomposition can be symbolically illustrated 
as 




+ 




(2.22) 



Hence, we can write the integral as a sum of box, left-vertex and right-vertex diagrams in 
the following way 

-^box + L.V. + R.V. = -^E.L.V. + -^E.R.V. > (2.23) 
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where 

-^E.L.V. 



da dj3 ^{ 



1 



1 



P j k^k + {x + 0P2]^ [k + pi + {x - 0P2 

which symboUcally means that 



2 ' 



(2.24) 



-^E.L.V. 




+ 



(2.25) 



while 



/E.R.V. = -|/ dadpd2k (-8)U 



k^{P + x-0\ 1 



+ -{l + a)iP + x-0 



P (x-0 j ^k+pi + {x-0P2]^ [k + {x-0P2 
which symboUcally means that 



, (2.26) 



Ie.r.v. = 




+ 



(2.27) 



We refer to Eq. (2.24) as effective left-vertex (E.L.V.) and Eq. (2.26) as effective right- vertex 
(E.R.V.) for simplicity and we consider the integrals separately. 

We start our analysis with E.L.V. and use Cauchy integration to integrate over a. A 
detailed analysis for the distribution of the poles is given in App. A.l. We are free to choose 

K 2 

to close either on the two poles corresponding to cutting the gluonic line, i.e. Cig = and 

I 2 

left fermionic line, i.e. a/ = ^^^jji^^g) ; or on the right and s— channel fermionic line, which 
we avoid. The integration over k is performed by using dimensional regularization. Then, 
the integral for E.L.V. reads 



-^E.L.v. = — 27ri / d^ I d^k ReSag + / dl3 d^k ReSaj , (2 
L-'O JO J-f-x Jo J 



.28) 



where 



Res 



1 



ag 



1 



-4 



(x - ix + 



+ 



W + x + iP + x + O 



2^ 



(x -I- i){x 




[x 




1 



\P{P + x-Ok 



2 ' 



2aP + x-0\k^ + W+--^0s 

(2.29) 
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and 



Vlpo — —A \ 



1 



1 



1 



-4- 



1 (/3 + X + 



1 

2^, 



1 



(2.30) 



with N = 2 — ejjv 
following relation 



2 + e/i^. In order to get the above expressions we have used the 



AB 



1 



1 

5 



(2.31) 



Using the fact that in dimensional regularization any integral without scale vanishes, the 
second line of Eq. (2.29) and the last term of Eq. (2.30) give zero. The ultraviolet divergence 
in k integral in both expressions is taken into account by renormalization and we are only 
interested in the finite part. 

Keeping all these remarks in mind, let us first calculate the gluonic pole contribution. 



-fE.L.V.,g = 




+ 



(2.32) 



The integration over k in the first term of Eq. (2.28) and using Eq. (2.29) gives 

27ri /-^-^ .J /3 1 , (/? + .T + (/3 + x + O 



-^E.L.V.,g = 4 



X 



dp 



+ 



(x + 0(x-0 2i{P + x-0 
sl3{l3 + X - ^y'"" 



xV{€uv) 



(2.33) 



Before integrating over f5, in order to simplify the calculation, one should stress that we are 



looking for the terms that contribute to the terms, i.e. the most singular terms. 

This corresponds to terms that are most singular at the limits of the integration. Hence, 
for the /E.L.v.,g integral we are interested in ^ terms which are singular at the lower limit, 
and pj^x-^ terms which are singular at the upper limit. In Eq. (2.33) there is no term 
that is proportional to ^. Thus, there is no contribution from ^ in /E.L.v.,g- The second 
singularity that should be considered is and for this type of singularity the integral 

reads 

/e.l.v., g = -4 -new) / dp{l5 + x- iY""-^ ■ (2-34) 



Integration over /3 is straightforward after this point. Hence, integrating over (3 and ex- 
panding the expression in the limit eiR^ (/y ^ 0, the finite part reads 



rfin. _ . 27ri 

-'E.L.V.,g — 



1 



?2! 



log'(e-x) 



(2.35) 



Let us emphasize that this contribution is coming only from the last term of the first line 
in Eq. (2.29) which is originated from the box diagram. 
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In the same way as in the calculation of the gluonic pole, the contribution from the 
fermionic pole, 



Ie.l.y.j 




+ 



(2.36) 



after k integration, gives 



-^E.L.v.J = 4 



27rz 



2e 



1 \V{^uv) 
(2.37) 

The most singular terms that we arc looking for in this integration are q , , , and -5-7^^ — 7 
terms which are singular at the lower and upper limit respectively. It is obvious from 
the expression that there are no pj^^j^^ type singularity. For ^j^^-e, ^^P^ singularity, the 
integral reads 

/E.L.V.,f =4— -r(e[;K) / dpiP + x-O'""-^- (2.38) 

^ S J -x-S, 

Again, integrating over j3 and expanding the result in the limit eir^ jjv ~^ Oj the finite part 
reads 

^RL.v.,f = 4^^1og2(20, (2.39) 

which is less singular than Eq. (2.35). Note that taken separately, each of these two 
diagrams which involve cutting the fermionic line lead to log^(^ — x)/{x — ^) terms, but 
this type of contributions add to zero at the end. 

A similar analysis can be made for the effective right vertex, E.R.V. Again two poles 
corresponding to cutting the gluonic line, i.e. oig = =^ and right fermionic line, i.e. oif = 

s(i3+x-^) ' considered. By using dimensional regularization the integral for E.R.V. is 
written as 



-^E.R.v. = — 27ri / dl5 I d^k Resag + / df5 djqk Resaf , 
Uo Jo J-f-x Jo J 



(2.40) 



where 



■{x-0^ 
1 



? x-^ 2e 



^2 _^ f^W+x+Qs 



+4- 



1 1 



1 



U2 ' 



(2.41) 



and 



. 1 rr 1 1 1 

/?ps = —4 < I 



(2.42) 
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with as above N = 2 — ejjv = 2 + eiR. 

Again, using the fact that any scaleless integral vanishes in dimensional rcgularization, 
one can immediately set the second line of Eq. (2.41) to zero. Moreover, one can see that 
the finite part of Eq. (2.42) vanishes totally since one of the terms is k independent and 
the second is scaleless. Hence, the only non- vanishing part of /erv after integrating over 
k reads 



-^E.R.v. = —4 



2m 



dp 



T{euv) 



<^IR 



(2.43) 

The most singular terms that may contribute to terms are ^ and terms 

none of which arc present in Eq. (2.43). The most singular term in the a; — >■ ^ limit is 
log ^) ^ j^Q]2cc /e.r.v. does not give contribution to '^'^ terms. 



x-i 



Thus, the above one-loop analysis shows that the only contribution to 



l°g'(g-^) 



terms 

come from Eq. (2.35) which is originated from the box diagram in the case of cutting the 
gluonic line around /3 + x — ^RiOin the phase space. We have thus shown that 



rdominant 
one loop 




27rz 1 \ 
-'^2!^°^(^-"^- 



(2.44) 



We would like to emphasize that the precision of our calculation does not permit us 
to fix the multiplicative coefficient a of (^ — x) under logarithm, i.e. Eq. (2.35) can be 
equivalently written as 



rdominant 
one loop 



, 27ri 1 , 2r //- M 

^ ^2! 



(2.45) 



The coefficient a is fixed to ^ by comparing the form of the log^(^ — x) terms in the exact 
one-loop result Eq. (1.7). Moreover, the shift ^ — > ^ — ze correctly takes into account the 
imaginary part of Eq. (1.7) leading to the following final formula 



rdominant 
one loop 



-4 



2m 



x-i + i€2\ 



log^ 



2e 



— le 



(2.46) 



This matching condition will be repeatedly used in calculations in higher loops and also in 
the resummation process. 

To conclude this section, we can thus state the first rule: 

(i) To extract the dominant behavior of the amplitude, it is sufficient to restrict ourselves 
to the contribution of the gluonic pole. 



This rule will be extended later after studying in detail the two-loop contributions. 
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3. One-loop in semi-eikonal approximation 



As explained in details in Sec. 2, the dominant contribution for x — >■ ^ is obtained from 
ladder-type diagrams at one-loop level. Thus, we now concentrate on the box diagram, see 
Fig. 6, and we show that the result that was obtained without making any approximation 
can be reproduced by using eikonal techniques applied to the left fermionic line of the box 
diagram. We do not rely on the Ward identities of the previous section, but rather expand 
the gluon propagator in the light-like gauge, using the fact that it will be considered to be 
on-shell. This will simplify much the analysis, similarly as in the case of small-x physics 
considered in Ref. [36]. The corresponding integral Ii reads 



Pi - 2CP2^ 



{X + C)P2 



Figure 6: The one-loop ladder diagram which contribute in the light-like gauge to the leading 

[as ln^(^ — x)]/ {x — ^) terms in the pcrturbative expansion of the DVCS amplitude. The p2 and _L 
momentum components are indicated. On the right, the dashed lines show the dominant momentum 
flows along the P2 direction. 



where 

(Num)i = tr{ hl,.[h + hMM + {x + hbujd'^" , (3-2) 

Ll=[ki + {x + 0P2]^ , S^=[ki+pi + ix-0P2]^ , Rl=[ki + ix-0P2]\ {3.S) 

and e = jll^^ +^i + {x- Ohha±- 

Wc now use an eikonal coupling for the left quark line, and we treat the gluon as 
soft with respect to this quark. Thus, in the incoming quark numerator ki + {x + C)P2 
is replaced by {x + C)P2- Furthermore, we treat this gluon as soft with respect to the 
s— channel fermionic line, thus working in the limit ai ^ 1. This leads to the fact that 9 
can be approximated asO = Note that this was checked in detail in Sec. 2, where we 

have seen that the dominant integration region corresponds indeed to the approximation 
ai < 1. 

Since the gluon is on mass shell, the dominant contribution from the gluon propagator, 
d^^ ^ when written in terms of gluon polarization vectors, is given by 

^"'^--i:^A)<A)- (3.4) 

A 

The numerator with the eikonal coupling to left fermion leads to the expression 

(Num)i = -2(x + 5^tr{ t>2l^Xh + - H i>i t>2 /(A)}(-ef;,)) , (3.5) 
A 
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which, after using the Sudakov decomposition of gluon polarization vector in pi gauge, 



-(A) 



a ^1 



(3.6) 



is rewritten as 



(Num)i = -2{x + E ( - ^27m[A + - M h h ^K-ef, 



Summing over the polarizations one gets 



(3.7) 



(3. 



Substituting Eq. (3.8) to Eq. (3.7), we arrive to the following expression for the numerator 



(Num)i 



-2{x + i) 
-2(x + 



tr<' h ( ^±1 - 2^ i*! ) [h + {x-0 h] h 



-2k j 



Pi 



1 + 



2(a;-g) ' 
/3i 



Pis 



2s 



(3.9) 



The appearance of 



, 2(x-g) 
Pi 



in Eq. (3.9) reflects the fact that the coupling of the right 

fermionic line is not the conventional eikonal coupling, since it takes into account some 
recoil effect. 

The denominators with on shell gluon (A;^ = 0) are 
Ll = ai{x + Os , Rl = -kl + ai{Pi+x-Os , S'^ = -g + {/3i + x - . (3.10) 

We use Cauchy integration to integrate over ai. The resulting expression for the residue 
at the pole ai = is 



Res, 



-4s 



1 + 



2(x-0 



Pi Ikiki-iPi + x- Os 

which leads to /i integral (see App. A.l for the limits of the /3i integral) 



4s- 



27ri 
X — ^ 



^—X POO 

dPi / dNki 
Jo 



1 + 



2(x-0 

Pi 



1 



1 



ki ki-iPi + x- 



(3.11) 



(3.12) 



Using the identity (2.31) for the last two propagators and taking into account the vanishing 
of the scaleless dimensionally regularized integrals, the expression (3.12) is 



L = -4 



in which also 



1 + 



2m 



/ dp^ / 
7o Jo 



dnki 



{p^ + x-Okl-{Pi+x-Os ' 



(3.13) 



is approximated by —1, in accordance with the ordering stated in 

Eq. (1.11) and illustrated by Fig. 6. This approximation is in accordance with the detailed 
analysis presented after Eq. (2.33). Integration over ki gives 

K-x 



h 



^ 27ri 
x-^ 



Jo 



dPiTieuv) 



1 



iPi+x-O 



{^-x-p^r'^ . 



(3.14) 
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From the regularized expression (3.14), one sees that in the /3i integration the dominant 
contribution comes from the region /3i around 0. Keeping this in mind, the finite part of 
the integral (3.14) can be determined from 



/r- = 4 



Iog(^ - a; - /3i) , 



(3.15) 



in which we already remove the regularization. 

Finally, after integrating over f3i and using the matching with the exact one-loop result, 
the dominant part of the one-loop diagram is 



rdominant 
one loop 



27ri 



1 



2C 



(3.16) 



in agreement with the result (2.46). 



4. Two-loop order 

Let us examine the next order in the perturbative expansion. There are many diagrams 
contributing but it can be shown that in the chosen gauge, the double box diagram domi- 
nates. 



4.1 Two- loop in semi-eikonal approximation 



Pi - 2^P2 



s 



%1 

R2 

Ri 

{X - P2 



V 



Figure 7: The two-loop ladder diagram which contribute in the light-like gauge to the leading 

[aj. ln^(^ — x)]/ {x — ^) terms in the perturbative expansion of the DVCS amplitude. The p2 and _L 
momentum components are indicated. On the right, the dashed line show the dominant momentum 
flows along the p2 direction. 

The analysis for one-loop case showed that the dominant contribution comes from the 
case where the gluon is on shell. So for the two-loop case we assume the same argument, 
i.e. both of the gluons are on shell. Moreover, assuming a strong ordering in \k^\ and Pi 
is natural, since the singularities to be extracted are leading double logarithmic ones. In 
practice, this means that we work in the approximation 



1^2! » 1^1 1 and x~^> l^il ~ > |a;-^-F/3i| ~ 1^2! 



(4.1) 
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which imphes, for on-sheh gluon (the fact that the gluon can be taken on-shell for the 
dominant contribution has been justified above), since we consider the gluons to be soft 
with respect to the s-channel fermion, that 



1 » \a2\ 3> |ai| 



(4.2) 



Note that the reverse ordering 1^2! <S |^i| would lead to a suppressed contribution, due to 
a non maximal number of coUinear singularities, which can be traced when evaluating the 
virtualities of the various loop momenta, as we will show in Sec. 4.2. 

As a consequence of these two assumptions we arrive to the fact that the coupling to 
the left fermionic line can be considered as eikonal coupling whereas the coupling to the 
right fermionic line is beyond eikonal approximation in a way that it takes into account 
the recoil effect (see Eq. (3.9) and the following remark.) Then, one can write the integral 
I2 for the two-loop ladder diagram as 

I J J daidpid2h J '^"2cZ^2C?2fc2(Num)2-^-^-^-^-^^^, (4.3) 
where the numerator is 

(Num)2 = -2(x + 0^ 
and the propagators are 



-2kj 
Pi 



1 + 



2(x-0 
Pi 



-2fe| 



P2 



1 + 



2(/3i+x-0 
P2 



2s 



(4.4) 



Ll 



ai{x + C)S , RI = -kj + ai{Pi + X - , 
a2{x + 0s , Rl = -kl + a2{Pi+p2 + x-0s, 
-kl + {Pi+P2 + x-0s . 



(4.5) 



We perform Cauchy integration over ai and a2 taking the residue at the pole /cf = A;| = 0. 
This leads to the result 



Res 



ai,a2 



4s- 



1 



X 



1 + 



2(x-0 
Pi 



1 + 



2{Pi +X-C) 
P2 



fcf {/3i+x- klkl -{P1+P2 + X- Os 



(4.6) 



The integral I2, Eq. (4.3), is written as 
h 



^ (27ri)2 f<-- 
4s- / dPi 

X-i Jq jo 



S,—x—l3i poo 

dP2 / li2^2 

Jo 

1 11 

X 



dok 



■2K1 



Pl+X-^klgg-{Pl+P2+X-0s ' ^^■'^^ 

with the limits of the P integration determined according to the discussion in App. A.l. 



Moreover, the expressions 



1 + 



Pi 



and 



1 + 



2{j3i+x-i) 
IS2 



are both approximated to -1 



according to the discussion after Eq. (2.33). Using the identity 



/*i2 r°° 

/ ^2^1 = / (^2^1 - / d2ki , 



(4.8) 
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which effectively shifts infrared divergences into ultraviolet divergences, one can use di- 
mensional regularization for the first integral in ki which then vanishes since it is scaleless. 
Thus, only the second integral contributes and the two-loop integral is written as 



-4s^ ^ / dpi / dP2 / 

Jo Jo 

1 11 



disfki 



i+x-^kikiki-i/3i+/32 + x- 



(4.9) 



Using the identity (2.31), 
1 1 



{f5i + P2 + x- Os 



(4.10) 



the first term vanishes in dimensional regularization since there is no scale. Then the 
integral I2 can be written as 



I2 = -4 



df3i / dP2 
Jo 

00 /"OO 

X / dNk2 



^1 + ,x - e /3i + /S2 + X - ^ 

/•OO poo 2 

/ dNk2 / div^i-^-n 

Jo Jkl '-^ '"^ 



kiki-{Pi + P2 + x-i)s 



(4.11) 



Integrating over k^ within dimensional regularization and taking only the finite term, we 
have 



rfin. 



^(2vri) 
x-i Jo 



dPi 



df32 



1 



1 



Pi+x-^l3i + P2 + x-^ 



/o 



kl-{/3i+ /32 + x- ^)s 



(4.12) 



Integration over ^2 can be performed in the same way. At this point we note that the same 
result can be obtained without invoking explicitly dimensional regularization, but using 
the method that is described in the App. A.2. Thus, using Eq. (A. 33) we get 



7-fin. 



(27ri 



X 



Jo 



1 



i-x j-i-x-fii 

1 



x^log2(^-x-^i-^2) . (4.13) 



It is straightforward to integrate over /3i and (52- Using the matching condition with the 
exact one-loop result, the integral I2 reads 



rfin. 
-'2 



-4 



{2m) 



1 



x-^ + ie4:\ 



log^ 



^ — X 



le 



(4.14) 



4.2 Detailed analysis of the suppressed diagrams at tw^o-loop 



In this section, we study in detail some of the two-loop diagrams in order to infer the 
minimal rules on which we will rely to then show that any diagram except the ladder-like 
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Pi - 2^P2' 




Figure 8: The two-loop subleading cross diagram. 



one are suppressed. These rules will be enough to also justify in the next section that only 
ladder-like diagrams contribute at any order. 

Let us first consider the cross diagram illustrated in Fig. 8. It reads 



1111111 



(4.15) 



' 7-2 7-2 C2 o2 d2 1,2 i.2 ' 
ij-^ 1^2 ^2 1 2 

with the numerator (Num) given by 

(Num) = tr|^2 A ( Y + ^ - ^^2^ ( Y + /3 + X - e);52^(x + Ohhi^ + Ohh } , (4.16) 
where the denominators are 

Lj = -kl + ai(/3i +x + Os, Ll = -fe + k2f + (ai + a2){Pi + ^2 + a; + C)s , 
^1 = -^2 + «2(/32 + x-i)s, Rl = -{ki + k2f + (ai + a2){Pi +P2 + x-^)s, 
= -ih + k2f + (1 + ai + a2)(/3i + ^2 + x-0s. (4.17) 

The ordering which leads to the dominant contribution is provided by a strong ordering 
both of transverse momenta and coUinear momenta, to extract the maximal logarithmic 
contributions, as 

1^2! > l^il and a; ~ ^ » > |/32| (4.18) 

(or 1^2! ^ l^il and a; ~ ^ ;» \P2\ ^ can easily check by inspection that any 

other ordering leads to less power of logarithms. Using the ordering (4.18), the residue is 



Res 



ai,a2 



-As 



1 + 



Pi 



1 + 



0: 



1 



P2 \Wi+x-0k2\x-0 
1 1 



fc2'(/3i +x-0k2^-{/3i + l32 + x- Os 

and the integral to be computed is 



(4.19) 



4s{2m) 



rt,—x—pi poo f^i2 

dPl / dP2 / ^2^2 / d2h 

JO JO JO 

111 1 



X-^k2'{x- k2' k2' -{Pl+(32+X- Os 



(4.20) 
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It is instructive to compare this expression witli Eq. (4.7). We explicitly see that integration 
over kl is different, since there is no A appearing which is the source of one power of 
log(^ — x). This is due to the fact that this cross diagram does not generate maximal 
collinear singularities. It thus shows that the net result can be neglected with respect to 
the dominant contribution ~ ^^^f-^^^^. 

The same reasoning applies to the ladder-like diagram which wc have discussed in 
Sec. 4.1. The Fig. 9 shows the virtualities of the various propagators in the two possible 
ordering in k^, keeping the usual ordering in /3j, namely x ^ \x — S,\ ^ 

\x — ^ + /3i \ ~ \P2\- The left diagram, where the ordering \k2\ 3> \ki \ is assumed, exhibits a 
maximal number of collinear singularities while the right one, with the opposite ordering 
1^2! ^ 1^1 1) does not. This justifies the kj^ ordering which was used in Sec. 4.1. 



kl + A 



kl 



1^ 



t hi 



kl 



kl + A 



kl 
kl 



f kl 
kl 



Figure 9: The two-loop ladder diagram which the two orderings. In both cases, the virtualities 
of the various propagators arc indicated. Left: natural ordering | fcj | ^ |^i 1 5 leading to the dominant 
contributions. Right: \ki\ ^ 1^2! leading to a suppressed contribution. 



Prom the above study, we can now infer the second guiding rule to extract the leading 
contribution in powers of log(^ — x), namely: 



(a) Each loop should involve a maximal number of collinear singularities, which manifest 
themselves as maximal powers of 1/^| for each i, after the ai integration according to 
rule (i). 



We now consider the diagram of Fig. 10, which involves the coupling of a gluon to 
the s— channel fermionic line. Closing the ccj contours on the gluonic poles, the fermionic 
propagators get virtualities whose order of magnitude are indicated on Fig. 10. Two limits 
are of interest in order to obtain the maximal powers of log(^ — x). The first one is the 
limit ^2 ^ that case, the number of collinear singularities originating from ki is 

too low, since there is a single propagator of virtuality fef (compensated by a similar fcf in 
the numerator), and this contribution is subleading. The second one is the limit fc| ^ M- 
In this case, the fact that the upper left fermionic propagator has a virtuality fc| + A 
where A = — (x — ^ + ^2)* lowers the level of singularity, again leading to a suppressed 
contribution. 

From this study, wc can now infer the third guiding rule to extract the leading contri- 
bution in powers of log(^ — x), namely: 
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M + A ki + A' 



k^ + A ki + A' 



kl l\ 



A 



k 



\ 



V 



Figure 10: A two-loop subleading diagram with a gluon attached to the s— channel fermionic line. 
The virtualities of fermionic lines are indicated. Left: ordering ^ ^ fc^. Right: ordering ]^ <C fci. 



(iii) Any coupling of a gluon to the s— channel fermionic line leads to a suppressed contri- 
bution. 

The last rule is obtained through the study of a diagram involving a fermion self- 
energy, of the type shown in Fig. 11. The key point here is to realise that the virtuality 
of the s— channel fermion is fcf + A, where /S. = —{x — ^ + I5\)s . The fact that A does not 
involve reduces the power of log(^ — x) when integrating over (32- 

From this study, we can now infer the fourth guiding rule to extract the leading con- 
tribution in powers of log(^ — x), namely: 

(iv) The diagram should he sufficiently non-local in order that the s— channel fermionic line 
involves the whole p2 flux. 



kl + A 



\k2 



kl " 



kl + A 
' ki 
kl 

kl 

ll 



Figure 11: A subleading two-loop diagram of abelian type. The virtualities of the propagators 
are indicated when closing on the gluonic poles. Left: ordering k^ 3> fci- Right: ordering k^ <S ^i- 



These four rules are sufficient to show that any non ladder-like diagram is suppressed, 
as we show now. The 3 diagrams of Fig. 12 are suppressed after applying rule (ii). The 5 
diagrams of Fig. 13 are suppressed after applying rule (iii). And the 4 diagrams of Fig. 14 
are suppressed after applying rule (iv). This last rule also excludes diagrams with virtual 
corrections on the gluon propagator. 
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Figure 12: Subleading two- loop diagrams violating rule (ii). 





Figure 13: Subleading two- loop diagrams violating rule (iii). 



Figure 14: Subleading two-loop diagrams violating rule (iv). 



5. All- loop analysis 



5.1 Beyond the two- loop order 

Based on the four rules formulated previously, it is now possible to justify that the con- 
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tributions to the maximal powers of 



only arise from the ladder-like diagram, at 



any order a" . For that, we arc using a recursive argument. At two loop, wc have seen 
that the diagrams with 3-gluon coupling are subdominant, since the powers of fcj are not 
maximal in that case. At three loop, the last missing building block, namely the 4-gluon 
vertex, appears. Since it is a contraction of two three-loop diagrams (which are already 
excluded) with one less propagator, this kind of vertex is also subleading. 

Thus, starting from the ladder- like diagram at order n — 1, let us dress it having in 
mind that we are looking for the maximal power of log(^ — x) , which will look ultimately 
like j^log2"(^-x). First, we are only allowed to consider abelian-like diagrams. Starting 
from a gluon which is attached somewhere on the right fermionic line, this line should end 
up on the left fermionic line: ending on the right would be too local (rule (iv)), and ending 
on the s— channel fermionic line would violate rule (iii). Finally, a crossing of any gluon 
line is not permitted since the rule (ii) would not be satisfied. Thus, we end up with the 
ladder-like diagram of order n. 



5.2 Computation of the n-loop ladder diagram 

The computation of the n-loop ladder diagram (see Fig. 2) is performed in full analogy 

with the one- and two- loop diagrams discussed in Sections 3 and 4.1. We again assume 
that all rung gluons in the ladder diagram are on the mass shell, and that there is a strong 
ordering in their transverse momenta and Sudakov variables a^, /3j, i.e. 



\kJ > \k 



> 



> 1^1 1 



1 » |q!,,| > la 



n-l 



» • • • > ai 



(5.1) 



X~^> \f5i\ r^\x-^\^\x-^ + f5i\^ |^2| >•••» \x-^ + /3i+/32 +/3n-l| ~ \Pn\- 

The above assumptions permit us to have eikonal coupling on the left fermionic line and on 
the right fermionic line the coupling goes beyond the eikonal coupling taking into account 
some recoil effects. Thus, the integral I„ for the n-loop ladder diagram is written as 



(S^\^ f f 111111 

-J j dax dPi d2^i ■ ■ ■ j dan d/3n d2^„ (Num)„-^ • • • JJ^^ 

where the numerator (Num)„ takes the form (compare with Eqs. (3.9) and (4.4)) 



(Num)„ = -2(x + 0' 



with the propagators 

Lf = ai{x + 
Ll = a2{x + S,)s 



-2k 



Pi 



1 + 



2(x-0 
Pi 

Pn 



-2^2 



1 + 



P2 

2{Pn-l + 



1 + 



2{Pi+x-0 

P2 

+ P1X-O 



Pn 



Rl 



-ki + ai{Pi + X - ^)s , 

-kl + a2{Pi+p2 + x-^)s, 



1 

(5.2) 



2s, (5.3) 



Ll = an{x + C)s 

c2 . 



-kl + aniPl -i \- Pn + X - , 



K + {Pl + ---+Pn + X-Os- 



(5.4) 
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Calculating the residue in aj, we get 

1 1 



1 



hi kiki-{ih + --- + Pn + x- Os 



(5.5) 



where each expression 



1 + 



^ ^ /32 



and 



1 + 



2{l3„-i+-+Pix-0 



IS 



approximated to -1 according to the discussion after Eq. (2.33). Then the integral /„ reads 

(27ri)" f^~^ fi-x-Pi Pn-i \ \ 



In = 4s(-l) 



On-l 



/ d/3,--- 
Jo Jo 

.oo 1 1 1 

/ d2k„--- ^^2^172 •••72 72 7^5— 

Jo Jo hi knkn- [Pl-^ 



/3l + X - e A + • • • + Pn-l + X-^ 



, (5.6) 



in which the limits of the /3j integrations are determined according to the discussion in 
App. A.l. Now we will use dimensional regularization to integrate over momenta k^. We 
use the identity (4.8) for the integral over ki. The first integral on the right hand side of 
the identity (4.8) vanishes. Using the same argument for n — 1 momenta and using the 
identity (2.31) for the product of the terms with momentum k^, the integral In can be 
written as 



-4^/ d^,... d/3n 
x-^ Jo Jo 

f^oo poo 







\ \ (5 7) 

M ]i^-l^n-{Pl + --- + Pn + X-Os' 



Integrating over fc^ , • • • , and only retaining the finite contributions from each integral, 
we have 



fin. 



n ri—x 



-4 



(27ri)' 
x-i Jo 



1 • • • 



i-x-Pi /3„_ 
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Jo 







00 -1 
,_.n-l ^ 

■n' 



X / d^u-iy 



Pl+X-^ Pi + ---+Pn+X-^ 
log-lfe2 



{n-iy.kl-{Pi + --- + Pn + x-Os' 



(5.8) 



The result of the last integral over k^ is given by Eq. (A.33) in App. A.2. In this way we 
obtain 



I; 



fill. 
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(2m 

x-i Jo 



/ dp,--- 
Jo Jo 



^-x-l3i /3„_ 



df3n 



1 



1 



/3i + .T - e Pl + ---+Pn+X-^ 



x(-i)"_log-(^-x-^i Pn). (5.9) 



Keeping in mind the remarks about the dominant region of /3i integrations after Eq. (3.14) 
and using the matching condition with the exact one- loop result, the Pi integrations lead 
to 



.fin. ^ . (2vrz)" 1 



log 



2n 



^ — X 



le 



(5.10) 
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5.3 The resummed formula 



The Eqs. (1.14) and (5.10) permits us to perform the resummation of the ladder diagrams 
and we obtain 



Kn=0 



X — ^ + i€ 



cosh 



£»log 



^ — X 



— le 



(x^-x). (5.11) 



The resummed to all orders formula Eq. (5.11) can now be included into the NLO coefficient 
function Eq. (1.4). The inclusion procedure is not unique and it is natural to propose two 
choices. The first case corresponds to modifying only the Born term and log^ part of 
Eq. (1.7) and keeping the rest of the terms unchanged. This corresponds to the following 
expression 

o2 



X — ^ + ie 



cosh 



L>log 



2^ 



le 



D2 

2 



9 + 3^-^ log 



— X 



(a; —x) . 



(5.12) 



In the second case the resummation effects are accounted for in a multiplicative way for 
Cq and Cl, i.e. the resummed formula takes the following form 

^2 



g\res2 



X — ^ + ie 



■ cosh 



I? log 



^ — X 



1-— <|9 + 3^^1og 



— X 



%)-{x^-x), (5.13) 



2tt ■ 



where D = 

These resummed formulas differ through logarithmic contributions which are beyond 
the precision of our study. 



6. Conclusions 

The resummation of soft-collinear gluon radiation effects allowed us to get a close all-order 
formula that modifies significantly the coefficient function in the specific region x near ib^. 
The measurement of the phenomenological impact of this procedure on the data analysis 
needs further analysis with the implementation of modeled generalized parton distributions 
and the discussion of specific observables. Let us just remind the reader that the region 
X = ±^ is crucial in the determination of beam spin asymmetries. 

We did not study the case of gluon GPD contributions to DVCS, which, although they 
are absent at Born order, are expected [27] to become important in the small ^ regime 
which will be accessible at high energies [37, 38] . 

Deeply virtual Compton scattering is but one of the exclusive processes giving access 
to GPDs. Our analysis could and should be applied to other processes too. The case of 
timelike Compton scattering is special since, thanks to the analyticity properties in Q'^, it 
has been shown [39] that a simple relation was relating its NLO correction to the one for 
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DVCS. The case of exclusive meson production is also very interesting, both theoretically 
and experimentally. The NLO analysis [23, 24] of the corresponding coefficient functions 
exhibit also a log^[(^ — x)/2^] behavior, both in the quark and in the gluon channels. It 
will be most interesting to see if our semi-eikonal analysis allows to resum these logarithms 
too. The quality of the present and near future data for vector mesons (and in particular 
p) electroproduction demands this analysis to be vigorously pursued. 

We did not study the effects of the running of a^. Also, a formulation of resummation 
in our exclusive case in terms of (conformal) moments is not yet available. This would 
generalize analogous resummation of inclusive DIS cross-section which were performed in 
terms of Mellin moments. We leave studies of these issues for a future work. 
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A. Appendices 

A.l Extracting the /Jj-ranges from the positions of the poles in a. 

In this appendix we show how the fii integration range is obtained from the study of the 
position of ai poles. 

Let us start with the qi pole. The tree denominators k\, Li and Ri possesses simple 
zeros in ai. Due to the ie prescription, their position in the complex-ai plane is governed 
by their imaginary part. For these three denominators, we respectively get 



kl: (A.l) 
Pi 

If 

Without loss of generality we can choose to close the ai contour in order to always avoid 
the pole due to Ri. Since the dominant contribution comes from closing on the pole of 
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, this gives the 4 following possibilities out of 6 
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(A.6) 

(A.7) 
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(A.9) 



Combining (A.4) and (A.9), the restriction due to the position of the 1/L^ pole in complex 
plane disappear. Similarly, one can combine (A.5) and (A.8). Thus one ends up with 

^ e{-Pi)e{Pi + x-o (A.io) 



Ri 
Ri 



-e{l3i)9i-Pi-x + 0, (A.ll) 



ki ) 

where we have used Eqs. (A.l, A.3). Finally, the Pi integral symbolicaly thus reads 

-2m / dPi Resai . (A.12) 





The same analysis can be applied for the 02 poles. Their position in the complex a2-plane 
is governed by 

kl : (A.13) 
and following the same line of thought we choose to close the 012 contour in order to always 
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avoid the pole due to R2- This gives the 4 following possibilities out of 6 
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(A.17) 
(A.18) 
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Combining (A.16) and (A.21), the restriction due to the position of the I/L2 pole in 
complex plane disappear. Similarly, one can combine (A.17) and (A.20). One thus ends 
up with 

^ e{-f32)e{l3i+(32 + x-0 (A.22) 



R2 
R2 



-e{l32)e{-l3i-(32-x + 0, (A.23) 



ki ) 

where we have used Eqs. (A. 13, A. 15). Finally, the ^2 integral symbolicaly thus reads 

-2m / dh Resaj • (A.24) 
Jo 

We apply the same analysis up to the q;„_i poles, which position in the complex Q;„_i-plane 
is governed by 

feLi ■■ -i^, (A.25) 

^"-1 • "TT^ — — (^-26) 

Pl + h Pn-l +X + i 

Rn-1 ■■ ^ '1 ^ 7- (A.27) 

/?i H h Pn-\ + a; - C 

and leads to the /3n-i integral of the form 

r^-x-Pi P„-2 

-1m / dPn-i Res„„_, . (A.28) 
^0 
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The last stage is achieved by considering the a„ poles. Since the additional fermionic prop- 
agator which joins the two photon vertices has a pole which position along the imaginary 
axis is the same as the one of their position is governed by 

kl: (A.29) 

Pn 

- ft + ...4+.-r 

and leads to the /3„ integral 



-27ri 



ri-x-pi p„_i 

/ dPnReSa„. (A.32) 

Jo 



A. 2 Some useful integrals 

In this appendix we show that 

y + A 



dy 



^ \ogP+^ A . (A.33) 



finite 



P+1 



Consider the integral 

Ip{A)= / -^dy, {AM) 
Jo 



which can be rewritten as 

s l ogP(At ) 
t + 1 

Its derivative reads 



Ipi^) = / 
Jo 



dt . (A.35) 



and thus, in the limit A — >■ oo , 

4(A)~-i^ + |Vi(A). (A.37) 
Now, based on the definition (A. 34), this derivative can be also expressed as 

from which one deduces that Ip{A) is UV finite. Thus, since from the definition (A. 34) it 
is clear that Ip is UV divergent, the divergency of Ip-i is exactly given by 

Vi(A)ldiv. = ^log^A> (A-39) 
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with no subleading logarithmic divergent terms, in order to compensate the first term of 
Eq. (A. 37). We can thus write symbolicahy 

and we now determine the dominant part of -^p(^) Ignite ' hmit A — >■ 0, which we 
denote as ip(A) . A direct calculation shows that 

7o(A) = log A - log A (A.41) 

so that 

7o(A) = -logA. (A.42) 
Thus, using Eq. (A. 37), this leads to 

7((A)~-i5|^, (A.43) 

and based on Eq. (A. 40) we obtain 

log^ A 

Ii(A) (A.44) 

/^(A)~-i^, (A.45) 
loe;^ A 

72(A) = -^|-. (A.46) 
By inspection, it is thus natural to make the ansatz 

Ip{A) = -CplogP+'A, (A.47) 
where cq = 1 . Combining Eqs. (A.37) and (A.47), we obtain the relation 

{p+l)cp=pcp-i, (A.48) 

and thus 

(A.49) 



The same analysis gives 



and thus 



which ends up the proof of Eq. (A. 33). 
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